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Abstract 

We develop the Bardeen-Cooper-Schrieffer (BCS) formalism for the superconductiv¬ 
ity of carbon nanotubes. It is found that the superconducting transition temperature 
T c of single-wall carbon nanotubes decreases exponentially with the increase of the 
tube diameter because the density of states near the Fermi energy is inverse propor¬ 
tional to the tube diameter. For the multi-wall carbon nanotubes, the Cooper paring 
hopping between layers enhances the superconducting correlation and increase the 
superconducting transition temperature, which is consistent with the experimental 
observation. 
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1 Introduction 


Carbon nanotubes (CNs) have promised a great potential application in nanotechnology. [T] 
In recent years, superconductivity has been observed experimentally in the 
single-wall carbon nanotubes (SWCN), the rope of SWCNs,|2j and the multi¬ 
wall carbon nanotubes (MWCN)[3] by the low-temperature transport mea¬ 
surements. For the individual SWCN and the rope of SWCN, the supercon¬ 
ducting transitions have been observed at a temperature below IK. [2] For 
the end-bonded MWCNs, the superconducting transition temperature reaches 
6 ~ 12iC,[3] which is much higher than that of the individual SWCN and the 
rope of SWCN. This implies that the interlayer coupling may enhance the 
superconducting correlation. [3] The investigation of the magnetic properties 
reveals that the superconductivity emerges at temperatures below 15 K for a 
small diameter tube (4 angstrom) embadded in a zeolite. [4] 
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From a theoretical point of view, it is still not very clear how to understand the 
mechanism of superconductivity in CNs. Generally, CNs may be regarded as 
a one-dimensional material (ID). The ID systems face some obstructions that 
prevent the emergence of superconductivity, such as Luttinger liquid states 
with a repulsive electron-electron interaction[5j and Peierls phase transition. 
However, one believes that the proximity effect could play an essential role in 
superconductivity state in CNs. [2] Another scenario on superconductivity of 
CNs is obtained by the effective low-energy theory for interacting electrons in 
metallic SWCN, which predicts that SWCNs can be described as a Luttinger 
liquid with an attractive electron-electron interaction. [BJ The breathing modes 
specific to CNs can be the origin of a strong electron-phonon coupling giv¬ 
ing rise to attractive electron-electron interactions. [6] A short-range attractive 
electronic interaction originates from the exchange of phonons, jj] The single- 
particle hopping between neighboring nanotubes in a rope is strongly sup¬ 
pressed because the different helical structure of nanotubes, in general, leads to 
the mismatching of the atoms between tubes. The intertube coherence is estab¬ 
lished mainly through the tunneling of Cooper pairs, which enhances the su¬ 
perconducting correlation. [7] Exactly speaking, CN is a quasi-one-dimensional 
system. The energy band contains some split bands near Fermi level. The inter¬ 
layer coupling in MWCNs could play an important role in superconductivity of 
CNs. Superconductivity in CNs should be intrinsic PIT] and the nature of su¬ 
perconductivity in CNs is still attributed to Bardeen-Cooper-Schrieffer(BCS) 
mechanism. [8] Thus, the BCS formalism for CNs is still an interesting issue. 

In this paper, we will develop a formalism of MWCNs based on BCS theory, 
and discuss the superconducting correlation and the critical temperature. In 
the section II, we will propose a BCS Hamiltonian of MWCNs in the tight- 
binding approximation. We will give a mean held solution of this Hamiltonian, 
the self-consistent equations on the superconducting energy gap and the su¬ 
perconducting transition temperature of MWCN in the section III. Finally, 
we will give a discussion and conclusion. 


2 BCS Hamiltonian 


Since the intrinsic superconductivity is observed in SWCNs one may believe 
that the effective electron-electron interaction is attractive to form the Cooper 
pairs in the SWCN. We consider a general BCS Hamiltonian of MWCN con¬ 
structed from a coupling multi-layer BCS Hamiltonian. [9] The single electron 
hopping between layers is estimated to be much weaker than the Cooper pair 
hopping between layers due to the mismatching of atoms between layers. m 
The Cooper pair hopping between layers dominates the interaction between 
layers. [2] Thus, the BCS Hamiltonian of MWCNs may be expressed in terms 
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of the coupling SWCN BCS Hamiltonian. 


H = Y,H i {k)+Y,H t j h 
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is the Hamiltonian of the lt\i tube. The Vg in the second term of Eq. (121) de¬ 
scribes the intratube effective electron-electron interaction. The Hijy in the 
second term of Eq. (P) is the coupling between layers, where the sum (£, £') 
runs only for the nearest layers in the tight-binding approximation. The layer- 
layer coupling Hamiltonian may be written as 

'y ) '^i,i , ^k,t,<r^'k,£,a^'k',i’,<r ( 3 ) 

k,k',a 


where the Vf ji in Eq.(j3]) represents the intertube electron-electron interac¬ 
tion between the £th and £'th layers, which induces the Cooper pairing hop¬ 
ping between layers and are measured relative to the corresponding screened 
Coulomb potential. The 4 <ff (cy i(r ) is the creation (annihilation) operator 
at the layer l with the spin cr, where the l = 1,2 ...M labels the layers. 
The Ei(k) = ee(k) — E F , where E F is the Fermi energy. We set E F = 0 
for convenience in the following. The q(/c) is the energy dispersion relation 
of the layer t for free electrons, which can be obtained by the tight-binding 
approximation {10] 
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and their wave vectors are 
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where a = 2.46A is the lattice constant of hexagons. The t is the intralayer 
electronic hopping constant. The SWCN is identified by the chiral vector 
C h = which length Lg = a^Jnf + mj + 7i £ rri£. The (ki,qg) in Eq.()5|) 

are quantum numbers of the £th tube. The sum J2k i n Eqs.([2]) and (J3J) runs 
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all states of the nanotubes i and t' in the Briullien zone, which includes the 
kg values within — jr < kg < jr and qg = 0, 1, ...Ng — 1, where Tg = \f3Lg/d rie 
is the length of the translational vector of the £th tube being parallel to the 
tube axis and Ng = 2 L‘j/ar , dn l is the number of hexagons per unit cell of the 
7 th tube. The dn t = gcd(2n + m, 2 m + n ) means the greatest common divisor 
of 2 ng + mg and 2 mg + ng. For comensurate MWCNs the atoms between lay¬ 
ers are matching such that all Tg are equal. [13J The energy band structure of 
SWNTs depends on the chirality of the tube. The SWCN may be metallic for 
ng — mg = 3 i or semiconducting for ng — mg ^ 3 i, where i is an integer. [TO] 


3 Mean field solution 


Based on the BCS mean field idea, factorizing the potential terms in Eq.((2]) 

and d3j) , Ck/,cr C k,t,a C k,£',a C k/',(T ~ { c k,l,<y C k,£,a) C k,i',^ C k,t',cr~\~C k ^ (y C k ^ & {Ck AA c k,l',a) ~ 
( C 1 e a c l i a)( c k,e',v c k,e',a), we can obtain the mean field form of Hamiltonian 

H = J2 ee(k)cl Aa c kAa - J2 a 2 c U,A/ a + h.c.) + ^ A \S ek (6) 
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where the pair correlation, Sg k , has been defined by 

Sgk \C k ,£,aCk,g,aT (7) 


where (. X ) is the grand canonical average of A". The superconducting energy 
gap may be expressed as [9] 

A^ = 22 VgSgk + 22 Vej'Se'ki (8) 
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to diagonalize the Hamiltonian in Eq.f|Hll. the Hamiltonian can be written as 
H B cs = E g + 22&{k){lk,f1k,t + (10) 





where the ground state energy is 


B 0 = EW*)-&(*) + j§y). HD 

and the quasiparticle energy spectrun, 

W*0 = ,/#) + A?. (12) 


where £ labels the tube, which may be regarded as different energy bands of 
quasiparticles. Making use of the Bogoliubov transform Eq. and considering 
the tight-binding approximation, the superconductng energy gap A^ satisfies 
a set of equations, which matrix form may be written as [9] 
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At zero temperature, the Eq. dlSl) may reduce to F(Ai) = ' .Ok 

In principle, numerically solving Eqs. (fl3]) and (1T5]) associated with Eqs.(ff|),([5|) 
and (TT21) . we can obtain the superconducting paring, superconducting transi¬ 
tion temperature and their relationships. The F(Ag) should satisfy Eqs. fflSl) 
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and (1151) . The nontrivial solution of A requires the determinant of the coeffi¬ 
cient matrix being zero. This way to solve numerically the BCS Hamiltonian 
is exact in the mean held level. 


3.1 Single-wall carbon nanotubes 


In order to give an analytic solution, let us consider the SWCNs first. We note 
that the density of states (DOS) of metallic tubes per atom can be expressed 
approximately near Fermi level as iromi 

P ( Ef ) - ± ( 1 6 ) 


where d is the diameter of the tube. Thus, the DOS near Fermi level depends 
on the diameter of the tubes. The equation of the superconducting energy gap 
may be written as 

1 = Vp(E F )F 0 ( A) (17) 


Since the Debye energy of the phonons hw B is estimated to be O.leV [7], we 
may neglect the effect of the energy band structure beyond the Fermi level. 
Thus, the function Fq(A) may be written as 


hw D 

Fo(A) = / 7FTW taah 


/Ve 2 + A 2 \ 
V 2 k B T ) 


(18) 


At the critical temperature, T = T c , the superconducting paring equals to 
zero, A = 0. The integration Eq. fl 18 D can be integrated 


F 0 ( 0 ) = In 


/ 2 e 7 ftujD \ 
\ * k B T c ) 


(19) 


where 7 ~ 0.5772 is Euler constant. Using Eq. (fT 6 ]h (TT7|) . and (TT9|l . we can 
obtain the superconducting transition temperature, 

, m 2 e 7 ( 7 r 2 td\ . . 

k B T c = - frw D exp — (20) 

7 T y 4aV J 

Interestingly, the critical temperature T c decreases exponentially with the in¬ 
crease of the diameter of the tube. We also obtain the superconducting paring 
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at zero temperature is A(0) = 2hu>r,exp (— Thus, the superconducting 
paring can be also expressed approximately as 

{ , A(0) 

A(0) - (2vrfc B TA(0)) 1 /2 e -^ for T « Tc 
3.06fc s T c (l - J) 1 / 2 for T^T C 

Since there is an energy gap at Fermi level for the semiconducting SWCNs, the 
superconducting state cannot occur in the semiconducting SWCNs. However, 
it has been found that the layer-layer coupling of MWCNs may induce the 
semiconductor-metal phase transition. jl2|[l4] Actually, most of MWCNs are 
metallic due to the interlayer coupling. H3H3 Thus, the superconducting state 
may occur in MWCNs. 


3.2 Double-wall carbon nanotubes 


For the double-wall carbon nanotubes (DWCN), the Cooper paring equation 
Eq. (USD can be reduced to [91TS] 

/ l-V lPl (E F )F 0 (A 1 ) -V lt2 p 2 (E F )F 0 ( A 2 ) \f A 1 
y — V2 ,iPi(Ef)F 0 (Ai) 1 — V 2 p 2 (Ep)Fo(A 2 ) J ^ A 2 

The effective electron-electron interaction may be assumed to be independent 
of the chirality of the tube, namely V\ = V 2 = V and let V 1:2 = V 2 ,\ = 
V±. At T — T c , the superconducting paring equals to zero, (A!,A 2 ) = 0, 
The nontrivial solution of Eq. (j22j) requires the determinant of the coefficient 
matrix being zero, in which substituting Eq. fllbp into Eq. (l22p . we can solve 
two solutions of Fq(0) of Eq. fl22p . Based on the hints of the experimental 
results [2|3] and the similar two-band BCS theory, [9] we may select the smallest 
solution of F 0 (0), 

= 8 -Vj)a |y(dl + dl) ~ A 2 (^ ~ <M 2 + WlEhl (23) 

Combining Eq. (fl9p . the superconducting transition temperature can be ob¬ 
tained 

k B T c =—hw D e- F ^\ (24) 

7T 

Thus, the interlayer hopping of the Cooper pair enhances the superconducting 
paring and increases the superconducting transition temperature, which is 
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consistent with the experimental observation [2lf3] . 


3.3 Multi-wall carbon nanotubes 

Similarly, we can generalize above way of the DWCN to the MWCN case. Sup¬ 
pose MWCN containing M layers, to facilitate the analytical solution of F 0 (0), 
we assume that V^p({E F ) = V F p F (E F ) = V(p(E F )) and Vij + ip£ + i(E F ) = 
Ve +1> tp e (E F ) = V±{p(E F )), where {p{E F )) = The (d) is the average of 

the diameter of MWCN, (d) = jjJ2edg. At T — T c , the determinant of the 
coefficient of the Eq. (1T3h can be reduced to the Chebyschev polynomial, which 
can be solved, j9] 


(25) 


(p(E F )){V + 214008(^)1’ 


where we has also selected the smallest solution of the Chebyschev polynomial 
like the DWCN case. [9] The superconducting transition temperature can be 
expressed as 



(26) 


It can be also seen that the Cooper pair tunnelling between layers will en- 
hacnce the superconducting correlation and increase the critical temperature. 
This agrees with the experimental observation. [2|3] This theoretical framework 
provides us an physical understanding of the superconductivity mechanism of 
CNs even though the result relies on the F 0 (0) solution selected based on the 
hints of the experiment result [2l3j and the multi-band BCS theory. [9j 

Experimentally, the transport measurement of the ropes of SWCNs indicates 
the superconducting state occuring in the temperature range of 0.1 ~ IK. [2] 
The investigation of the magnetic properties for the small radius SWCN (4 A) 
reveals that the T c can be estimated to be lower than 15A'.[3J The transport 
measurements of the Au/MWCN/Au junctions indicate that the T c reaches 
to 6 K ~ 12K. [3] The theoretical study of the zigzag and armchair SWCNs 
by the Millan’s formula and the VASP software package gives that the T c is 
very small, 10 -9 /l for the armchair tube (5,5), but 55 ~ 75 K for the zigzag 
tube (5. 0).[T5] 


As an example of this formalism we list several typical results of the supercon¬ 
ducting transition temperature in the table I and compare the first-principle 
calculation[3] and the experimental results. [21415] where we use the parameters 












Table 1 

The superconducting transition temperature of several CNs. 


CNs 

(n,m)@(n,m)... 

UK) 1 

UK) 2 * 

UK)' 

SWCN 

(Tl) 

0.729 

57 

< 15 

SWCN 

(5,5) 

9.23 x 10' 4 * * * * 

3 x 10" 9 

0.1 ~ 1(SWCN rope) 

SWCN 

(10,10) 

6.46 x 10” 10 

2.6 x 10 _3 ° 


DWCN 

(5,0)@(10,0) 

2.24 



DWCN 

(5,5)@(8,8) 

0.043 



TWCN 

(5,0)@(10,0)@(15,0) 

0.29 


6 ~ 12(MWCN junction) 

TWCN 

(5,5)@(8,8)@(12,12) 

0.0069 




1 our results; 2 the first-priniciple calculation[15j: 3 experimental results. [2l3l4| 


t = 2.5eV, V = 1.2eV, and V± = 0.8eV, Hlud = 0. leTT. [7j For SWCN the first- 

principle calculation predicts the (5,0) tube being metallic due to the tt* — a 
coupling. |l5j However, in the tt electronic tight-binding approximation the 
(5,0) tube is semiconducting.[TO] Thus, we chose a small diameter tube (4,1) 
(~ 4 A) for comparision with the results of the experimental observation [1] and 
the first-principle calculation. [15] The result we obtain is T c = 0.729 K, which 

agrees qualitatively with the experimental result. [2ll3ll4] For SWCN (5,5), we 
estimate the T c = 9.23 x 10 -4 /i, which is the same order to another theoretical 
results, T c = 1.5 x 10 _4 /l .[16] In their theory the electron-phonon interaction 
is taken into account in detail within the BCS framework. For MWCN the 
Cooper pair tunnelling between layers also enhances the Cooper pair corre¬ 
lation and increase the superconducting transition temperature. This is also 
consistent with the experimental observation. m 


4 Discussion and conclusion 


Physically, the intratube electronic hopping constant t, the effective electron- 

electron interaction Vi and the intertube electron-electron interaction 

could depend slightly on the chirality and the layer number of the tube. PUS] 

Thus, We cannot obtain exactly an analytical solution of T c . Nevertheless, 

above formalism Eqs.([l])~ f|T3|) still provides a guideline for the numerical 
study. From the theoretical point of views, we give the BCS formalism on 

MWCN and the formula in Eq. (Effij) tell us that the Cooper pair tunnelling be¬ 
tween layers enhances the superconducting correlation and T c , which is cconsis- 
tent with the experimental results [2ll3] even though we make some approxi¬ 
mations in Eq. fl26p for obtaining an analytic formula. Actually, for MWCN 
the effect of the layer number on T c competes with the effect of the diameter 
of MWCN, which can be seen from Eq. (12611 . However, the effect of diarne- 
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ter dominates T c for the MWCNs that have the layer number more than 4. 
Hence, the small-diameter double-wall or triple-wall carbon nanotubes could 
have higher T c . 


On the other hand, the study of the thermodynamic variables starts from 
the partition function of the system, Z = Tr(e -/3// ), where (3 is the inverse 
temperature. The free energy may be expressed in terms of F — — f3 In Z. 
Using Eq. flTOjh we can obtain 


F 



ln(l + e -fo( k ))- ee (k)+Z t (k) 


A 2 1 
2 &(*)_ 


(27) 


The other thermodynamic variables can be obtained by the thermodynamic 
relationships between the free energy and the thermodynamic variables, such 
as the electronic entropy, S — ^§f), and the specific heat C — T .[8] 

In summary, we have developed a BCS formalism of superconductivity for 
MWCNs, which give that the superconducting transition temperature T c of 
SWCNs decreases exponentially the tube diameter. The T c is very small (less 
than IK for the diameter less than lnm). The Cooper pair tunnelling between 
layers enhances the superconducting correlation and increases the superson- 
conducting transition temperature. The interlayer coupling increases the num¬ 
ber of the transverse channels, which suppresses ID characteristic of CNs due 
to screening the repulsive electron-electron interaction. Hence, the supercon¬ 
ducting long-range order emerges favorably in the small-diameter double-wall 
or triple-wall carbon nanotubes . 
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